Theory of Josephson Arrays in a Resonant Cavity 



E. Almaafl 

Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556 

D. Strou43 

Department of Physics, The Ohio State University, Columbus, Ohio 43210 
(Dated: February 2, 2008) 

We review our previous work on the dynamics of one- and two-dimensional arrays of underdamped 
Josephson junctions placed in a single-mode resonant cavity. Starting from a well-defined model 
Hamiltonian, which includes the effects of driving current and dissipative coupling to a heat bath, we 
write down the Ifeisenberg equations of motion for the variables of the Josephson junction and the 
cavity mode. In the limit of many photons, these equations reduce to coupled ordinary differential 
equations, which can be solved numerically. We present a review of some characteristic numerical 
results, which show many features similar to experiment. These include self-induced resonant steps 
(SIRS's) at voltages V = nhQ/{2e), where 57 is the cavity frequency, and n is generally an integer; 
a threshold number Nc of active rows of junctions above which the array is coherent; and a time- 
averaged cavity energy which is quadratic in the number of active junctions, when the array is 
above threshold. When the array is biased on a SIRS, then, for given junction parameters, the 
power radiated into the array varies as the square of the number of active junctions, consistent with 
expectations for coherent radiation. For a given step, a two-dimensional array radiates much more 
energy into the cavity than does a one-dimensional array. Finally, in two dimensions, we find a 
strong polarization effect: if the cavity mode is polarized perpendicular to the direction of current 
injection in a square array, then it does not couple to the array and no power is radiated into the 
cavity. In the presence of an applied magnetic field, however, a mode with this polarization would 
couple to an applied current. We speculate that this effect might thus produce SIRS's which would 
be absent with no applied magnetic field. 

PACS numbers: 05.45.Xt, 79.50.-|-r, 05.45.-a, 74.40.-|-k 

I. INTRODUCTION 

The properties of arrays of Josephson junctions have been of great interest for nearly twenty yearsi These arrays 
are excellent model systems in which to study phenomena such as phase transitions and quantum coherence in two 
dimensions. For example, if only the Josephson coupling energy is considered, and if the self-inductance and mutual 
inductance of the array plaquettes are neglected, the Hamiltonian of a two-dimensional (2D) array of Josephson 
junctions is formally identical to that of a 2D XY model [see, e. g. Ref. In the presence of such inductive 
effects, this XY description needs to be modified, and several generalizations, which include these effects, have been 
proposedi^i^*^ Arrays sometimes appear to mimic behavior seen in nominally homogeneous materials, e. g. high- Tc 
superconductors, which often behave as if they are composed of distinct superconducting regions linked together by 
Josephson coupling^ Finally, the arrays are of potentially practical interest: they may be useful, for example, as 
sources of coherent microwave radiation if the individual junctions can be caused to oscillate in phase in a stable 
manner. 

Recently, our ability to achieve this kind of stable oscillation, and coherent microwave radiation, was significantly 
advanced by a series of experiments by Barbara and coUaboratorsi^iSiifliiLiS These workers placed two-dimensional 
underdamped Josephson arrays in a resonant microwave cavity. The presence of the cavity caused the junctions 
to couple together with remarkable efficiency, resulting in a highly efficient conversion of the injected d. c. power 
into a. c. radiation. Even more surprising, this efficiency is achieved in underdamped arrays, which according to 
conventional wisdom, should be especially difficult to synchronize, since each such junction is both bistable and 
hysteretic. Numerous workers have attempted to explain these experiments li^tiiii^ii^ 

In this brief review, we summarize our own model to describe this kind of phase locking^Siil In Section II, we 
present the model Hamiltonian. In Section III, we describe the equations of motion resulting from this model, in the 
"classical limit," as defined further below. In Section IV, we give a brief description of the relevant numerical results. 
Finally, in Section V, we present a concluding discussion. 



2 



II. MODEL HAMILTONIAN 



We consider a 2D array of TV x Af superconducting grains placed in a resonant cavity, which we assume supports 
only a single photon mode of frequency fl. The array thus has {N — 1){M — 1) square plaquettes. There are a total 
of X Ny horizontal junctions, where Nr^ = N — 1 and Ny = M. A current / is fed into each of the M grains on 
the left edge of the array, and extracted from each of the M grains on the right edge. Thus, the current is injected in 
the X direction, with no external current injected in the y direction. A sketch of this geometry is shown in Fig.^ We 
also introduce the terminology that a "row" of junctions, in this configuration, refers to a group of Ny junctions, all 
with left-hand end having the same x coordinate, and all being parallel to the bias current. One such row is indicated 
by the dashed lines in Fig. ^ 

We write the equations of motion for the grain variables (phases and charges). We express our Hamiltonian in the 
formsiSiii 

H = Hphoton + Hj + He + Hcurr + Hdiss- (1) 

Here Hphoton is the energy of the cavity mode, expressed as 

Hphoton = Kl ^a^a + , (2) 

where and a as the usual photon creation and annihilation operators. Hj is the Josephson coupling energy, which 
takes the form 

Hj^-Y^E-lcos{^,,), (3) 

{ij) 

where E^^ is the Josephson energy of the (ij)*^ junction, and 7,,, is the gauge-invariant phase difference across that 
junction. E'-^ is related to If^ , the critical current of the (ij)*'' junction, by Ef^ = hlfj/q, where q = 2|e| is the Cooper 
pair charge. He is the capacitive energy of the array, which we write as 

^'^ = ^E9'(^^"'k'^^«.' (4) 

where is the inverse capacitance matrix, is the number of Cooper pairs on the i*'' grain, and g = 2e is the 
charge of a Cooper pair (we take e > 0) . We write ^ij as 

7y ^(t)^- <j)j - [(27r)/$o] / A • ds = ~ - A,j, (5) 

J ij 

where is the gauge-dependent phase of the superconducting order parameter on grain i, $0 = hc/{2e) is the flux 
quantum, and A is the vector potential, which in Gaussian units takes the formMiiS 



A(x, t) = Vihc^)/m («W + «^ W) E(x), (6) 

where E(x) is a vector proportional to the local electric field of the mode, normalized such that Jy d'^x|E(x)p = 1. 
Again, f2 is the resonant frequency of the cavity mode, and V is the cavity volume. The line integral is taken across 
the (ij)*'' junction. The phase factor Ay is then 

Aij = gij{a + flt), (7) 

where 



/fic2 (27r)3 



9v-\l—^ / E-ds (8) 



characterizes the effective coupling between the (ij)*'' junction and the cavity. 
One can also define a frustration for the /^*'' plaquette by the relation 



2tt 

plaquette 
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where the sum runs over the bonds in the /i plaquette. For the present case, 

/M=/r"' = ^(«+«^) E 9^,■ (10) 

plaquette 

If there were an appHed magnetic field normal to the array in addition to the cavity electric field, then the frustration 
would have an additional contribution 

jrnag ^ ^j^^^ 

where $^ is the magnetic flux through the plaquettei^S 

Note that gij is very sensitive to the experimental geometry. For example, if the cavity has the form of a par- 
allelepiped with edges L^, Ly, and L^, where > ^ Ly, then lowest mode in this cavity is a TE mode with 

frequency — nc^l/Ll + l/Ly] the corresponding value of is given 

2 2 32c s . . 

where is the cosine of the angle between the E of the resonant mode and the vector ds) . In the geometry of Ref. 
1^ Q/{27r) w 100 GHz, a cavity of this shape gives rise to gij 0.001,^^ though it is obviously very sensitive to both 
array and cavity geometry. 

The driving current and dissipation may be incorporated as followsiiS*ii The driving current is included via a 
"washboard potential," Hcurr, of the form 



f^jext 



q ^-^ 

(y>ll* 



where / is the driving current injected in the x direction into each grain on the left edge and extracted from the 
right edge, the sum running over bonds in the x direction. Dissipation is included by coupling each 7^ to a separate 
collection of harmonic oscillators with a suitable spectral densityi2i*22i2i24 

B^ss^Y.^^'^ (14) 

where the sum runs over distinct bonds (zj), and 

{Pa,ij ) 



jjdiss \ 



fa 



1 ifai-Y 



(15) 



The variables Ua^ij and Pa.ij of the a oscillator in the (ij)* junction, are canonically conjugate, and ma,ij and Wa,y 
are the oscillator mass and frequency. If the spectral density, Jij{uj) oc |a;|, the dissipation is ohmic.^--^- We write 
such a linear spectral density as 

Jiji^) = TT" 1^1 Q('^c - ^), (16) 

ZTT 

where ujc is a high-frequency cutoff, Q{uJc~uj) is a step function, and the dimensionless constant — Ro/Rij, where 
Rq = h/ (4e^) and i?y is a constant with dimensions of resistance, which turns out to be the effective shunt resistance. 

III. EQUATIONS OF MOTION 

The equations of motion are obtained by introducing operators a = an + iaj and = an — iaj . These satisfy the 
commutation relation [an, aj] = i/2, which follows from [a, a^] = 1. In terms of these variables, Hphoton — fin{a^ + aj) 
and 7y = 0; - (f>j - 2gijaR. 
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The time-dependence of the various operators appearing in the Hamiltonian is now obtained from the Heisenberg 
equations of motion. These are readily derived from the commutation relations for the various operators in the 
Hamiltonian Q. Besides the relations already given, the only non-zero commutators are [nj,(j)k] — —iSjk and 
[Pa,ijTUi3,ke] ~ ^oi,i3 Sij^ki, where the last delta function vanishes unless (ij) and {k£) refer to the same junction. 

Using all these relations, we find, after a little algebra, the following equations of motion for the operators (pi, rii, 
ail, and a/: 

i 1 



l-J2{C-%n„ (17) 

3 

^■^ X! ^^^(^i -(t^i- '^guan) 
I 

Text 1 



(18) 



Ua.iifa.il H — (^(f)^ - (f>i - 2giiaii) 

I a ^ ' ' 'ma,il{U}a,ilj 

ciR = fl ai, (19) 

J7]J_ jext ^ 

ai = -il qr + ^9v^ sm{(j)i ~ - 2gijaR) ^ gy 

+ T.irT. I /".^^ ""-^ + ^^"^'""r (<^^ - - 25.,«fl) 1 ■ (20) 

Here, the index I ranges over the nearest-neighbor grains of i. In writing these equations, we have assumed that the 
only external currents /f^* are those along the left and right edges of the array, where they are [cf. Fig. 

Eqs. p7|l - H2U() are equations of motion for the operators aR, ai, rij, and (or 7j). In order to make these equations 
amenable to computation, we will regard these operators as c-numbers. This should be reasonable when there are 
many photons in the cavityii^^ 

The equations of motion for the harmonic oscillator variables can also be written out explicitly, but are of no 
interest; we instead eliminate those variables and incorporate a dissipative term directly into the equations of motion 
for the other variables. Such a replacement is possible provided that the spectral density of each junction is linear 
in frequency, as above. In that case^iSiSiiSSiSiSi the oscillator variables can be integrated out. The effect of carrying 

out this procedure is that one should make the replacement (yfa,ij Ua.ij + lij^ ~^ ^F"^ T'y wherever 

this sum appears in the equations of motion. Making this replacement in the equations of motion and simplifying, 
we obtain the equations of motion for Uj and a/ with damping: 

ai = -n aR + ^g,j^ sin(7y ) ^ 5»j + XI ^^^^'^'i' ^^^^^ 

Once, again, the index j is summed only over the nearest-neighbor grains of i. Equations (|17(l . p9|l . H21|l. and H22|l 
form a closed set of equations which can be solved for the time-dependent functions 7^, rii, o-r and a/, given the 
external current and the other parameters. 

To express these equations of motion in terms of suitable scaled variables, we introduce a dimensionless time 
r = tqRI'^/h = ijjrt, where R and /"^ are suitable averages over Rij and Ify We also define the other scaled 

variables % = VL ^ I ^ j^, V, ^ orj = ^2'r:^aRj, = \[^^9tj, C'ij = i^rRCiy The last 

relation involves the capacitance matrix dj. We assume that this takes the for m^^i^^ dj — {Cd + Zi Cc)5ij — 
Cc (<5ij+x -I- i5ij_x -I- ^ij+y -I- <5ij-y), i. e., that there is a non-vanishing capacitance only between neighboring grains 
and between a grain and ground. Here Zi{— 4) is the number of nearest neighbors of grain i, Cd and Cc are respectively 
the diagonal (self) and nearest-neighbor capacitances, and x and y are unit vectors in the x and y directions. The 
corresponding Stewart-McCumber parameters are (3c — lUt-RCc and Pd = LOrRCd- 

In the above equations, we have introduced Vi = qJ2ji^~^)ij''^j' which is the potential on site i. The integral of 
the electric field across junction (ij) is Vij = Vi — Vj — 2gijD,ai. 

Carrying out these variable changes, and after some algebra, we find the dimensionless equations of motion 

^^ = V,, (23) 




1 . . . ^ 

- ^/ - 2gjiaH) + -~-{Vi -Vi - 2gui}ai) 
Rii I 



(24) 



Via J 



(25) 



and 



ai = —Clan + 



9rj 



I^j sin( 



Hij 



{'y>ll* 



(26) 



where the dot refers to a derivative with respect to r. These equations do not include the self-magnetic fields produced 
by the currentsi^i^i^ However, they can be generalized to include external currents with both x and y components, 
and non-square primitive cells. Also, note that the frustration parameter defined earlier is now time-dependent. 



IV. SOME NUMERICAL RESULTS 



We now describe some numerical results obtained by solving Eqs. (|23|l - (|26f) numerically,^- using the adaptive 
Bulrisch-Stoer methodSl as further described in Ref. 1161 . For simplicity the coupling constants gij were assumed to 
have only two possible values, gx and gy, corresponding to junctions in the x and y direction respectively. This 
assumption should be reasonable if there is little disorder and the resonant mode has long wavelength compared to 
the array dimensions. 

The calculated IV characteristics for the case g^ ^ 0, jjy — 0, with driving current parallel to the x axis are shown 
in Fig. 121 The array is taken to have 10 x 4 grains, with capacitances /3c — 20 and /3d — 0.05 (independent of junction), 
gx = 0.012, and ft = 0.41. The critical current through the (ij)*'' junction is Ifj = 1 + Ay where the disorder Ay 
is randomly selected with uniform probability from [—A, A], and A = 0.05. The product IfjRij is assumed to be 
the same for all junctions, in accordance with the Ambegaokar-Baratoff expressioni^^ The calculated IV's are shown 
as a series of points. The arrow directions of the arrows indicate whether the curves were obtained under increasing 
or decreasing current drive, or both. The horizontal dashed curves are voltages where self-induced resonant steps 
(SIRS's) are expected, namely {V)r /{NRIc) = tt, where {V)t is the time-averaged voltage (dotted lines are guides 
to the eye). Each nearly horizontal series of points denotes a calculated IV characteristic for a different number of 
active rows Na, and represents Na x Ny (horizontal) junctions sitting on the first integer {n = 1) SIRS. The calculated 
voltages agree well with expected values given by dashed horizontal lines. The long straight diagonal line segment 
represents the ohmic part of the IV characteristic with all rows active. (The corresponding segments for other choices 
of Na < 10 are not shown). Besides the integer SIRS's, there are a few fractional SIRS's, similar to what is seen for 
Shapiro steps in conventional underdamped junctions>22, 

Fig. O shows the IV characteristics for three different arrays, each with all rows in the active state: (i) a 40 x 1 (full 
curve), (ii) a 40 x 2 (dotted curve) and (iii) a 40 x 3 (long-dashed curve). Each array has the parameters gx = 0.015, 

= 0.49, Pc = 20, Pd — 0.05 and A = 0.05. Once again, the arrows denote the directions of current sweep. The 
horizontal dot-dashed curve shows the expected position of the SIRS corresponding to Na = 40 [V/ {NxRIc) = 
The curves show that all three arrays have qualitatively similar behavior. First, if the array is started from a random 
initial phase configuration, such that / = I/Ic > 1 -I- A, and / is decreased, then all the rows lock on to the Na = 40 
SIRS. Secondly, if / is further decreased, the Na — 40 active state eventually becomes unstable and all the junctions 
go into their superconducting states. Finally, if / is increased starting from a state in which the array is on the 
Na = 40 SIRS, the SIRS remains stable until I reaches the critical current for the various rows, and the IV curve 
becomes ohmic. 

In Fig. 0] we plot the time-averaged energy E{Na) = {a\ + aj)r in the cavity for three different arrays: 40 x 1 
(stars), 40 X 2 (circles), and 40 x 3 (squares). In all cases, / = 0.58, and the other parameters are the same as those 
of Fig. O Below a threshold value of Na, (which we denote Nc and which depends on Ny), the active rows are in 
the McCumber state (not on the SIRS's). In this case, E{Na) is small and shows no obvious functional dependence 
on Na (see inset). By contrast, above threshold, E{Na) is much larger and increases as iV^. Fig.^shows that, when 
Ny is increased at fixed gx, Nc decreases. Precisely this same trend is observed when we increase gx while holding 
Ny fixed (and was observed in our previous ID calculations with increasing gx). Thus, the relevant parameter in 
understanding the threshold behavior appears to be NyCjx- 
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For the 2D arrays, one can introduce a define a Kuramoto order parameter {rh)r for the horizontal bonds by 



where Na is the number of active rows, Ny is the number of horizontal junctions in a single row and the sum runs 
over all the active, horizontal junctions. For the parameters shown in Fig.^ it is found, as in ID,^^ that {rh)T ~ 1 
for Na > Nc while {rh)T ^ 1 for Na < Nc- This behavior occurs because, for this choice of parameters, none of the 
active junctions arc on a SIRS when Na < Nc- 

For the case — 0, g^j ^ 0, i'n our geometry (with square primitive cells), we have not been able to find any value 
for (jy for which a SIRS develops. This behavior is easily understood. In this geometry, with current applied in the x 
direction, there is little power dissipated in the vertical junctions and no resonance is induced in the cavity. 

It is no surprise that the cavity interacts only very weakly with the vertical junctions. From previous studies of both 
underdamped and overdamped disordered Josephson arrays in a rectangular geometry (see, e. g., Refs. Isol and lSlT) . 
it is known that when current is applied in the x direction, the y junctions remain superconducting, with {V)r ~ 0, 
while the x junctions comprising an active row are almost perfectly synchronized, with (rx) ~ 1. 

If there were an an external magnetic field perpendicular to the array, we believe that SIRS's would be generated for 
gy 7^ 0, even if g^ = 0. In this case, as mentioned earlier, there would be a non-zero magnetic-field-induced frustration 
jmag jgq^ As a result of this magnetic-field-induced frustration, there would be nonzero voltages across the y 

junctions, as well as supercurrents in these junctions. 

It is of interest to compare these 2D results explicitly with those of ID arrays. In Fig. [S] the IV characteristics 
of a 10 X 1 array having coupling constant 5a;;ioxi = 0.0259 with those of a 10 x 10 array with coupling constant 
5a;;ioxio — 0.00259. The other parameters are the same for the two arrays: (l — 0.41, /3c — 20, f3d = 0.05 and 
A = 0.05. The expected positions of the SIRS's [at V/{NRIc) = il] are indicated by dashed horizontal lines. Indeed, 
the two sets of IV characteristics are very similar. This indicates that the crucial parameter is Nyg^, not Ny and ijx 
independently. The slight extra flatness in the 10 x 10 IV's probably occurs because the individual junction couplings 
in the 10 x 1 array are 10 times larger than those in the 10 x 10 array. Similarly, the differences in the "retrapping 
current" in the two sets of curves (i. e. the current values below which the McCumber curve becomes unstable), are 
due to the fact that, for a given value of A, the 2D arrays are effectively less disordered than the ID arrays, since the 
average critical current for a single row has a smaller rms spread than the critical current of a single junction in a ID 
array. (For other discussions of the effects of disorder, see, e. g., Refs. [32 33 34). 

In both the 10 x 10 and the 10 x 1 array of Fig.|31 the width of the SIRS's varies similarly (and non-monotonically) 
with the number of active rows. This behavior distinguishes our predictions from those of some other models, "'^'^'^4 
where a monotonic dependence of SIRS width on Na is found^^ and the cavity is modeled as an RLC oscillator 
connected in parallel to the entire array. 

In Fig. El the reduced time-averaged cavity energy E = (a^ + a'j)r is plotted vs. I = I/Ic for both arrays of Fig. 
|31 under conditions such that all rows are active. This plot is obtained by following the decreasing current branch. 
When the 10 x 10 array (with g(ioxio) — 0-1 5(ioxi)) locks on to the SIRS, E jumps to a value which is approximately 
two orders of magnitude larger than in the 10 x 1 array, even though Nyg^ is the same for both arrays. The reason 
is that though the width of the steps is controlled primarily by Nyg^, the energy in the cavity is determined by the 
square of the number of radiating junctions, which is is 100 times larger for the 2D than the ID array. 



The equations of motion described lead to a transition from incoherence to coherence, as a function of the number 
of active rows Na- There is also a striking effect of polarization: the transition to coherence occurs only when the 
cavity mode is polarized so that its electric field has a component parallel to the direction of the current flow. 

The numerical results closely resemble the experimental behavior, showing the following experimental features: 
(i) self-induced resonant steps (SIRS's) in the IV characteristics; (ii) a transition from incoherence to coherence above 
a threshold number of active junctions; and (iii) a total energy in the cavity which varies quadratically with the 
number of active junctions when those junctions are locked onto SIRS's. But there may be some differences, as 
discussed further in Ref. 0. 

The 2D theory bears many similarities to the ID case and elucidates why the ID model works so well. These 
similarities occur because, in a square array, only junctions which are parallel to the applied current couple to the 
cavity. Thus, as in ID, the 2D model leads to clearly deflned SIRS's with voltages proportional to the cavity resonant 
frequency. However, there are some numerical results specific to 2D. For example, whenever one junction in a given 




(27) 



V. DISCUSSION AND SUMMARY 
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row is biased on a SIRS, all the junctions in that row phase-lock onto that same SIRS. In addition, when the array 
is biased on a SIRS, E(Na) is much larger in 2D than in ID, for the same value of the coupling parameter cjxNy. 

A key difference between ID and 2D is the effect of polarization: when the cavity mode is polarized perpendicular to 
the applied current, it does not affect the array IV characteristics. We believe that, if the array were frustrated, e. g. 
by an external magnetic field normal to the plane of the array, there would be a coupling even when the cavity mode 
is polarized perpendicular to the current. It would be of interest to carry out calculations for the model described 
here, to confirm this effect. 

We conclude with a brief discussion of ways in which our model might be generalized to include some effects which 
are presently not taken into account. First, of course, one should include the fact that the cavity, as well as the 
junctions, has a finite damping (finite Q). Second, the real cavity fields are not spatially uniform within the array, 
as assumed in the calculations presented here. We speculate that the effect of these two generalizations would be 
to inhibit the array synchronization. Third, all real cavities have more than one resonant mode; these other modes 
may be relevant in some experimental circumstances. Fourth, we have treated the operators as c-numbers, i. e., have 
neglected quantum effects related to their non-commutativity. These quantum effects will certainly be relevant in 
some circumstances. Finally, the present treatment of damping within individual Josephson junctions is also carried 
out in the classical limit (resistively and capacitively shunted Josephson junction); in the quantum regime where the 
individual junction variables need to be treated quantum-mechanically, this treatment will need to be modified. 

Most of these generalizations can probably be carried out in a straightforward manner. Thus, we believe that the 
present model has most of the essential physics underlying the SIRS's seen in experiments. Of the omitted effects, we 
believe that the quantum effects, if included, are most likely to produce qualitative changes, because they could lead 
to intriguing entanglement between quantum states of the array and of the cavity^ 
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Horizontal junction 




Vertical junction 



FIG. 1: Sketch of the array geometry in our model. There are (M x N) superconducting islands (black squares), making 
[{M — 1) X + (A'^ — 1) X M] Josephson junctions (crosses). An external current J*^^* is injected into each junction at one end of 
the array and extracted from each junction at the other end. The array is placed in an electromagnetic cavity which supports 
a single resonant photon mode of frequency O. The dashes denote a "row" of junctions, which is perpendicular to the current 
bias and is comprised of horizontal junctions. 
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FIG. 2: Calculated currcnt-voltago characteristic for a 10 x 4 array with cavity frequency Q, = 0.41, capacitance parameters 
/3c = 20 and Pd = 0.05, disorder parameter A = 0.05 and junction-cavity coupling in the horizontal direction gx = 0.012. 
The horizontal dashed lines show expected voltages for the various SIRS's. These correspond to different numbers of rows of 
horizontal junctions in the active state. Arrows denote that the given IV was taken for increasing or decreasing current. 
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Current l/l^ 

FIG. 3: Calculated current- voltage characteristics for a 40 x 1 (full line), a 40 x 2 (dotted line) and a 40 x 3 (long-dashed lino) 
array, all with = 0.015, tl = 0.49, /3c = 20, l3d = 0.05 and A — 0.05. The horizontal dot-dashed line shows the expected 
SIRS position. Note that as the array width increases, the smallest I at which all the active junctions phase-lock on the SIRS 
also increases. Hence, increasing the array width and increasing have a similax effect. 
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FIG. 4: Time-averaged scaled energy E in the resonant cavity as function of active number of rows for a 40 x 1 (asterisks), a 
40 X 2 (circles) and a 40 x 3 (squares) array with driving current 7 = 0.58. All the other parameters are the same as those of 
Fig- El Inset: an enlargement of the IV characteristics near the synchronization threshold, on a logarithmic vertical scale. The 
threshold number of active junctions for synchronization decreases with increasing array width. 
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FIG. 5: IV characteristics for a 10 x 1 array (*) and a 10 x 10 array (o). The 10 x 1) array has parameters Qxaoxi = 0.0259, 
O = 0.41, /3c = 20, [3d ~ 0.05 and A = 0.05. The expected SIRS positions of the SIRS's are marked by horizontal dashed lines. 
The 10 X 10 array has 3a;,ioxio = 0.00259, and the other parameters are the same as for the 10 x 1 array. The IV characteristics 
are shown for both increasing and decreasing current drive. 
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FIG. 6: Time-averaged reduced cavity energy E, for a 10 x 1 array and a 10 x 10 array for the same choice of array parameters 
as in Fig. |S] The calculations are carried out on the decreasing current branch with all rows active, for the 10 x 10 array is 
10 times smaller than that of the 10 x 1 array. 



